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Chiral perturbation theory for K 
on the lattice 



7r + 7r° decay in the continuum and 



Maarten Golterman a and Ka Chun Leung 3, * 

a Department of Physics, Washington University, St. Louis, MO 63130, USA. 

We use one-loop chiral perturbation theory (ChPT) to compare lattice results for the K + — > 7r + 7r° decay 
amplitude with the experimental value. Three systematic effects: quenching, finite-volume effects, and the use of 
unphysical values of quark masses and pion external momenta can be investigated. We find that the corrections 
help in explaining the discrepancy between lattice and experimental results. We also discuss the relation to Bk ■ 



1. Basic theory 

In the Standard Model, the decay K + — ► 
7r + 7T° is induced by the four-fermion operator 

Shl^d-L dhl^dh- In ChPT, O4 is represented by 



4 = a 27 r«(E^St)^(S^Et)/ 



(1) 



at 0(p 2 ) with the tensor rjy projecting out the 
AS = 1, AI = 3/2 components, and a large 
number of 0(p 4 ) operators, each with its asso- 
ciated coefficient, constructed from the unitary 
E- field and the quark- mass matrix M |Q . A sim- 
ilar representation can be constructed for O' = 

SlJ^cIl SLl^dL, which is related to K —K° mix- 
ing. At 0(p 2 ) this amounts to just changing the 
tensor r]y . O4 and O 1 are components of the same 
27-plet under SUl(3), and therefore the same co- 
efficient a 27 appears in both O4 and O' . 

The matrix element (tt + tt \O4\K + ) to one loop 
can be obtained by calculating loop diagrams 
from the 0(p 2 ) operator, which gives rise to chi- 
ral logarithms with known coefficients, and which 
depend on a cutoff A, and tree-level contribu- 
tions from 0(p ) operators. Since we do not 
have enough information to determine the 0(p 4 )- 
operator coefficients, we will estimate the size of 
one-loop corrections by setting all these coeffi- 
cients to zero and choosing A to be 770 MeV or 
1 GeV, which is generally believed to be the en- 
ergy scale below which physical effects of more 
massive hadrons can be absorbed into the 0(p i ) 
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coefficients. The sensitivity of physical quantities 
to these different values of the cutoff is taken as an 
estimate of the systematic error associated with 
the lack of knowledge of the 0(p A ) coefficients. 

2. Continuum result 

The one-loop calculation of (7r + 7r° |C>4 1^ -1- ) 
in ChPT (for m, = 0) was first undertaken 
in Numerically, the real- world result for 

(7T + 7r°|04|A" + ), with m u = m c i ^ m s , and ra v = 
136 MeV, m K = 496 MeV and /„ = 132 MeV, is 



(7r + 7r |O 4 |A' + ) = 



12ia 27 



(m 2 K - ml] 



1 



V2f* 

0.63, A = 1 GeV\ 
0.36, A = 770 MeV J ' 

We see that the one-loop correction is fairly large. 

3. Lattice results 

It is possible to extract a related matrix ele- 
ment on a lattice with spatial volume 1? from a 
computation of the correlation function 



(3) 



c(t 2 ,h) = <o|^+(; 2 v (t 2 )o 4 (ti)A-(o)|o) 

t2»ii»0 c -g 2x (t 2 -t 1 ) c -m Jf t 1 x 

(0|7r+7r |7r+7r°)(7r+7r |C>4|A: + )(A+|A-|0) 
(7r+7r°|7r+7r )(A-+|A-+) ' 

where 7r + (t) = 7r + (ir, i), ir + = 7r+(0), etc. 
and E^tt is the energy of a state with two pi- 
ons at rest in a finite volume. (One also needs 



2 



kaon and two-pion correlation functions.) Ex- 
cept for rriK = 2m v , an "unphysical" matrix ele- 
ment is obtained since all external mesons are at 
rest. Also, in current lattice computations un- 
physical masses (degenerate quark masses which 
are heavier than real- world ones), the quenched 
approximation and, of course, finite volume are 
used. We will use quenched ChPT H to cal- 
culate (TT+ir \O 4 \K + }i uenched to one loop. The 
result can then be compared to eq. @ to get 
an estimate of all these systematic effects. It 
should be noted that the coefficients a 27 and 
a 27 , the tree-level meson decay constants / and 
f q (q for quenched), etc., are in principle not 
equal. As for finite volume, we have: 1. the 
difference between finite-volume and continuum 
values for operator coefficients are exponentially 
small in L (ESL) ||; 2. spatial momentum in- 
tegrals J (%J? S f{k 2 ,m T! ) are replaced by discrete 
sums -p- J2k=^-fi f(k 2 , m>n),n € % 3 (for periodic 
boundary conditions). If the integrand is regu- 
lar, the sums are equal to the continuum inte- 
grals with corrections ESL. Otherwise, there are 
additional power corrections in L^ 1 H: 



L 3 ^ k 2 



d 3 k /(fc 2 ,m^) 
(2tt) 3 k 2 



0-22578 Id/ 
— /(0> ™*) - Jz (0, m * 



(4) 



ESL. 



We will not consider ESL corrections. We cal- 
culated the unphysical amplitude for degenerate 
quark masses, in which case O4 and O' do not 
couple to the ?7'-meson, and hence there are no 
contributions from the "//-parameters" 6 and a 

1- 

Details of the calculation can be found 
in g]. The diagram in which the two pi- 
ons produced from the if + -decay strongly 
rescatter leads to power-like finite-volume ef- 
fects. Suppose the strong-interaction vertex 
acts at t s with t% < t s < ti. The essen- 
tial part of the expression for this diagram is 



L~ 3 EE Jfas exp [-2(y/k* + ml - m n )(t s - t x )]. 

For k = 0, this gives [to, — ti)/L 3 which can be 
resummed into the tree-level result, and thus 
produces = 2m^ + I/(2L 3 / 2 ) in the ex- 



ponent in eq. (|). For k f_ 0, we get terms 



oc l/(y m 2 + k 2 — m n ) = ( y m 2 +k 2 + ra^)jk 2 . 
These give rise to sums like eq. (Q) , which results 
in L^ 1 and L~ 3 corrections. We will discard any 
excited-state contributions ||. 

The collected one-loop corrections for the un- 
physical matrix element (7r + 7r°|O4(0)|.K" + ) are 



2 / 2 

Tfi I HI 



(5) 



relative to the tree-level value 2Aia 2 7m^L 3 / V%f 3 
for the full theory, and 



? / 9 

777 / 777 



(6) 



relative to the tree-level value 2Afia\ 7 m 2 x L 3 / \/2f 3 
for the quenched theory. 

F{m K L) = 17.827/(m 7r L) + 127r 2 /("i 7r L) 3 (7) 
is the finite- volume correction. 

4. Numerical examples 

To one loop in ChPT the physical matrix ele- 
ment and the unphysical one from quenched lat- 
tice computations (after extrapolation to the con- 
tinuum limit) are related by 

u 
f 



{■K + AOM\K + ) phvs =Y^[ J -j 



rrr K - m; 



2M 2 



(7r+7r |O 4 (0)|^ 



, quenched 
I unphys ' 



with 



Y = 



0.089, A=l GcV 
-0.015, A=770 MoV 



f 



Ml 



o 1 M i 
-31og^ 



F(M„L) 



(8) 



(9) 



The conversion factor Y embodies all the one- 
loop corrections. The a 27 and / ratios in the 
prefactor are not known. We will arbitrarily set 
them equal to one. It remains one of the un- 
certainties that cannot be resolved within ChPT. 
Mk and refer to values computed on the lat- 
tice while rriK and m w refer to real- world values. 
When we apply the formula to the lattice data 
of (in which the mass-squared-ratio prefactor 



3 



was already taken into account), we get values 
which are shown in fig. 1, along with the origi- 
nal data (for which the error bars are statistical 
only). The error bars on our points come from 
varying A and A q independently, and do not con- 
tain the statistical errors. We have eliminated 
points with smaller physical volume or at which 
M„ > 770 MeV m m p . At all points, Y < 1 
and each "corrected" amplitude lies below the 
corresponding original one. The one-loop results 
reduce the discrepancy between the lattice data 
and the experimental result. However, one- loop 
effects are rather substantial, and two-loop cor- 
rections can probably not be neglected. For more 
discussion of all uncertainties involved in these 
estimates, we refer to |H|. 
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a 2 7 and / drop out, unlike in eq. (g). At M n = 0.4 
GeV, M^L = 8 and F n = 132 MeV, X ps 0.6. 

5. Relation to Bk 

One-loop results for Bk with m u = ^ m s 
can be found in ||. The ?/-double pole con- 
tributes and hence the result depends on <5 and 
a for m u = rrid 7^ m s . 

The ratio 



Tl = Sk (n + n°\0 4 \K+)/ (K°\O'\K ) 



(12) 



in which a 27 drops out is of some interest. In 
the physical case, the one-loop correction is 60% 
to 95% of the tree-level value in magnitude, de- 
pending on the values of the cutoff. This calls 
into question the reliability of ChPT in this case. 
In the unphysical case, 



srjfull rj-f quenched 

unphys unphys 



(13) 



3i 
71 



1 + 



Ml 



(4ttF w )2 



3 log 



Ml 
A 2 



F(M„L] 



The typical magnitude of the one-loop corrections 
for the data of § is 10-15%. 

We would like to thank Claude Bernard and 
Steve Sharpe for discussions. This work was sup- 
ported in part by the DOE. 
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Figure 1. Open symbols: data from || (squares: 
16 3 x 25 (or x33), (3 = 5.7; octagons: 24 3 x 40, 
(3 = 6); crosses: including the correction factor 
Y. The constant c = 2v / 2/(G^ sin# c cos^c). 



A similar relation involving instead the unphys- 
ical unquenched matrix element is 



(^ + 7r |O 4 
with 

X = 



Iphys 



X 



2M2 

full 



(n + n»\O,(0)\K + Y u :; hys , (10) 
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F{M„L) 
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